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We study experimentally a coiling instability of cylindrical multilamellar stacks of phospholipid 
membranes, induced by polymers with hydrophobic anchors grafted along their hydrophilic back- 
bone. Our system is unique in that coils form in the absence of both twist and adhesion. We interpret 
our experimental results in terms of a model in which local membrane curvature and polymer con- 
centration are coupled. The model predicts the occurrence of maximally tight coils above a threshold 
polymer occupancy. A proper comparison between the model and experiment involved imaging of 
projections from simulated coiled tubes with maximal curvature and complicated torsions. 
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I. INTRODUCTION 



Coiling is a common occurrence in vastly different systems and on many scales, ranging from carbon nano-tubes 
|jj and DNA molecules to telephone cords and tendrils of climbing plants ||. The framework most commonly 
used to model coiling in these and other systems is that of an elastic rod. Buckling of the central line is explained by 
showing that converting twist to writhe lowers the elastic energy H. 

We study the coiling of cylindrical stacks of lipid membranes H], called myelin figures, interacting with an am- 
phiphilic polymer. This polymer has several hydrophobic side groups attached along a flexible hydrophilic backbone, 
which it inserts into the membranes in order to shield them from the surrounding water j(| . The membranes compos- 
ing a myelin figure are in a two-dimensional liquid state, and therefore cannot support twist. Application of torque 
on the cylinder simply leads to flow of material around the tube. Thus the interplay between twist and writhe cannot 
explain the coiling observed in our experiments. 

The existence and coiling of myelin figures has been observed as far back as 1854 ]7j] ( see also ||). More recently, 
similar shapes have been observed during the hydration of a surfactant by brine PJlCj]. Coiling of myelin figures 
has also been reported for a system of egg- yolk lecithin (ll] [ll| . This system is a mixture containing many different 
lipids, having a variety of tail lengths and degrees of saturation Jl4| . In another experiment it was shown that a 
binary mixture of dimyristoylphosphatidylcholine and cardiolipin forms single and double helices in the presence of 
calcium ions p5[ . In these two studies, it was claimed that: (i) the energy gained by surface adhesion overcomes the 
energy cost of bending a tube, and (ii) the tighter the coil, the longer the line of contact between tubes becomes, thus 
increasing the area of contact. 

In contrast, our experiments clearly show that surface adhesion is negligible in our system. In order to account 
for coiling in our experiments, we present a simple model, in which we assume that polymer molecules locally induce 
spontaneous curvature. The coiling instability results from a coupling between local polymer concentration and 
membrane curvature. Within such a framework, hollow tubes may also undergo other shape instabilities such as 
pearling ||,[l6 17 1. However, the constraints imposed by the geometry of the cylindrical membrane stack prevent them 
from occurring. Our model predicts the occurrence of maximally tight coils above a threshold polymer concentration. 
Indeed, only maximally curved coils were formed in our experiments. A theoretical analysis shows that virtual slices 
through maximally curved helices can be very similar to the observed images. 

Section |l| describes the materials and methods used in our experiments. Section III describes the types of coils 
observed, and shows that adhesion does not play a major role in these experiments. Section IV presents the simple 
model that explains the major findings. First, a heuristic argument is given as to why such an approach works. Then 
the full mean-field calculation is given. Finally, section ^ outlines the difficulties in comparing quantitatively the 
model and the experiment. 
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II. MATERIALS AND METHODS 



Vesicles were made of stearoyl-oleoyl-phosphatidylcholine (SOPC) with Cig alkyl chains. The polymer used is 
hydrophilic dextran (MW 162,000 g/mol) functionalized both with palmitoyl alkyl chains and dodecanoic NBD chains 
as fluorescent markers. The hydrophobic anchors, distributed statistically along the backbone (about 1 alkyl chain 
per 25 glucose units) are C\§ long. On average there are about 4 persistence lengths between consecutive anchors. 
Therefore, the extension of each polymer molecule on the two-dimensional membrane is much larger than its extension 
into the third dimension. Events were observed by phase contrast microscopy and recorded on video. For fluorescence 
imaging the NBD markers were excited with argon laser illumination, and observed with a cooled CCD camera. 

Samples were prepared by drying a 0.5-1.0/il droplet of SOPC dissolved in a 4:1 chloroform-methanol solution 
(7.35mg/ml) on a glass slide. The sample was then closed from the top and sides, and hydration was effected by 
injecting a polymer solution of known concentration, c p , into the closed cell. Some time after hydration, a variety 
of self-assembled structures formed, including myelin figures. All the structures observed in these experiments were 
still connected to a lipid reservoir, allowing exchange of material, including both lipids and polymer molecules. 
Experiments were conducted at room temperature, well above the solid-liquid transition for SOPC. This allowed free 
diffusion of anchored polymers along the membranes. 

We stress that while the polymer concentration in solution, c p , is known, we do not control the surface concentration 
on the bilayers. The slow evolution of some of the structures we observe is consistent with a possible variation of this 
concentration over time. 



III. RESULTS 



Electron micrographs of cross sections of myelin figures reveal that they are rodlike lyotropic liquid-crystalline 
structures containing a large number (hundreds to thousands) of concentric cylindrical membranes separated by thin 
hydration layers The smectic order in these stacks of membranes is not ideal, as many defects are present. The 
outer radius of a myelin hgure can reach tens of microns, while the radius of the water-core may be of the order of, or 
smaller than optical resolution (0.2/i). The myelin figures, which are connected at one end to a large lipid reservoir, 
are continuously elongating. The rate of elongation of a myelin figure lies in the range of — O.Sfim/ sec. Throughout 
the experiment, polymer molecules continue to anchor from the surrounding solution. Hence we assume that there is 
continuous accumulation of both lipid and polymer molecules. 

Hydration of a patch of lipids by a polymer solution of small c p results in the formation of myelin figures, which 
display a clear tendency to straighten over lengths many times larger than their diameter (see Fig. Q). As c p is 
increased, myelin figures become more floppy and curved. For large enough values of c p , a writhing instability sets in 
and tubes bend, forming irregular structures (Fig. ||), single coils (Fig. ||a) and double helices (Fig. The type of 
coiled structure formed, depends in great part on the dynamics of the formation process. By far the most common 
event is for the tip to begin to curve in upon itself forming a seemingly irregular ball-like structure (Fig. ^|). However, 
when the ball is large enough, some sort of ordering can be seen (Fig. ||e). Coils that start forming at their bases 
usually evolve into nearly ideal single helices. The uncoiled part of the tube leading to the tip is pulled in and wound 
around to form the next loop |l8|. In cases where the coiling starts in the middle, it begins as a hairpin that rotates 
around itself, and the instability then proceeds in both directions (Fig. ||). This may result either in the formation 
of a tightly packed sphere or a double helix. At the site where the instability nucleates, the tube goes from nearly 
straight to maximally curved. The instability then propagates from this site to the rest of the tube (Fig. § and ||). 
The observed evolution may be due to a gradual change in the concentration of polymer molecules on the membrane. 

All the coiled structures we observe are maximally curved already as they form and do not tighten up gradually, 
unlike the experiments reported by Sakurai et al. |ll|]l2f| . In quantitative terms, this means that the curvature of the 
tube central line, C, is C « l/^o, where rg is the radius of the tube. 

In previously reported cases coiling of myelin figures was attributed to surface adhesion. This is clearly not 

the case in our system, as demonstrated by the experiment illustrated in Fig. |^. The tip of a myelin figure, in the 
process of coiling, adhered to an air bubble, which we then moved. Movement of the bubble stretched the coil (Fig. 

-c), until the latter reached a configuration in which all self-contact was lost (Fig. ||c). Note also that the coil is 
stretched more or less homogeneously. Upon detachment from the bubble, the coil retracted as if it was an ordinary 
spring (Fig. |d-g). 

Had surface adhesion been the dominant mechanism, one would not expect the response of a coil to mechanical 
stretching to be homogeneous, but rather to come apart at the site of weakest contact. Furthermore, adhesion cannot 
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create a restoring force. Thus, if the coil is stretched open so that no contact sites are left, the tube should "forget" 
that it was coiled. If the force exerted on the end of the tube is then released, adhesion would induce coil formation 
starting at one point and propagating to the rest of the tube (similarly to the original formation process) . This is in 
stark contrast with our experimental observations. 

Another piece of evidence against adhesion in our system is provided by the presence of many other structures in 
our experiments that come into contact with one another, but do not adhere (see e.g. Fig. ^|). 

In order to understand better the role that the polymer plays in the coiling phenomenon, it is important to know 
its location in the membrane stack. For this purpose, fluorescence imaging was used. The results are shown in Fig. [?]. 
As can be seen, the fluorescence intensity through the slice is concave. This is what one expects for a homogeneous 
polymer distribution throughout the stack, since in that case, the intensity should be roughly proportional to the 
thickness of the tube in the microscope slice. From this we infer that the polymer is present, in significant quantity, 
throughout the myelin figure. Images of hollow tubes, where the polymer sits predominantly on the outer layers, have 
a convex fluorescence intensity profile through their cross sections. 

IV. THEORETICAL MODEL 

We now present a simple theoretical model that captures most of the key experimental observations. For simplicity, 
we regard the system as if it were in equilibrium, effectively ignoring the slow evolution of the observed structures. 

Consider a stack of concentric cylindrical sheets. We represent each bilayer as two square lattices (in the spirit 
of lattice-gas models), corresponding to the outer and inner monolayers. Each site represents a patch of membrane 
having area a 2 , approximately the size of a polymer molecule performing a two-dimensional random walk on the 
membrane. Each lattice site has two degrees of freedom associated with it: the local mean curvature, H, and a 
binary occupation variable, <r, that takes the values one or zero when the site is occupied or unoccupied by a polymer 
molecule, respectively. To estimate the area of a site we assume that the hydrophilic backbone is in a good solvent 
(water) under semi-dilute conditions. This gives a radius of gyration, a w 40 — 80 nm. 

The energy of the system is a sum of the curvature energies of the sites: 2kH 2 for a vacant site, and 2k! (H — Hq) 2 
for a site occupied by a polymer, n and k' are the local bending rigidities of a monolayer without and with an attached 
polymer, respectively, and Hq is the spontaneous curvature induced by the polymer. We assume that Hq > 0, since 
the addition of polymer tends to bend membranes into shapes with higher curvatures. For the purpose of this simple 
model, the exact molecular mechanism responsible for this is not important. Possible mechanisms include the entropic 
pressure of the polymer backbone, or the incommensurability between the anchors and lipids. By convention, the 
curvatures of the inner and outer monolayers of the bilayer have opposite signs at the same position. 

A crucial assumption is that polymer molecules can diffuse along the membrane, since the membrane is in a fluid 
state. In addition, we assume that the fluctuations of each membrane in a stack are severely restricted by the presence 
of its neighbors. As a result, there is a strong correlation between the transversal fluctuations throughout the tube. 
The myelin figure can thus be regarded as a flexible rod, having a circular cross section everywhere along its axis. 
The experimental pictures indeed display unchanging circular cross sections within experimental error. 

Based on these assumptions, we developed a model that predicts that a high enough polymer concentration on 
the membrane can shift the equilibrium state, from a straight tube to a maximally curved one. We first present a 
heuristic argument to show that if the spontaneous curvature, Hq, is large enough, the free energy of a bent tube 
may be lower than that of a straight one. This approach may give a more intuitive understanding than the detailed 
calculations that follow. 

A. Heuristic Arguments 

In order to find the equilibrium state of a tube, we have to calculate its free energy. This free energy depends on 
the curvature of its central line, C, and on the polymer concentration. 

Consider one cylindrical bilayer of length I and circular cross section of radius r, with the same average polymer 
concentration, p$, on both sides, po is defined as the number of polymer molecules on a monolayer, n, divided by 
the total number of sites on a monolayer, i.e. po = na 2 /A, where A is the total area of the membrane segment. Let 
us calculate the free energy cost of bending the bilayer into a portion of a coil with central line curvature C in three 
steps. First, we bend the tube while keeping the distribution of polymer around the tube homogeneous. Next, we 
allow the polymer to diffuse from regions of lower curvature to regions of higher curvature (Fig. ^|), and finally, we 
consider the entropy of mixing. 
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The energy of a bent cylindrical bilayer with a homogeneous polymer distribution is Eh om = ^tom ^hom > wnere 
^hom an< ^ ^hom are ^ ne ener gi es of the outer and inner monolayers, respectively. According to our model 

£w(C) = 2/W J dA[H±(C,<l>) - H*]* + 2(l-po)« / dA[i/±(C»] 2 , (1) 

where is the angle around the tube, H ± = ±5 + 1 ^ c °^^ ^ is the local membrane curvature, and dA = 

dfa[l + Cr cos fa). For simplicity, we shall assume here that k' — n. In the full model we allow the possibility that the 
presence of the polymer affects the local bending rigidity (i.e. k' k). 

For our geometry the total mean curvature obeys / dAH^iC^fa — ±2ttI, independent of C. Thus the cost of 
bending the cylindrical membrane, keeping the polymer distribution homogeneous is: 

^Et m = Et m {C) - Et m (0) = 2nJdA [H±(C, fa? - H±(0, 0) 2 ] . (2) 

From this it is obvious that there is always an energy cost to bend the tube in this way. However, the price is 
independent of the value of Hq. 

Next, we take into account inhomogeneities in the distribution of polymer around the tube. Such inhomogencitics 
reduce the energy, if polymers move to regions of membrane curvature closer to Hq in both the outer and inner 
monolayers. Rewriting Eq. (|l|) with p being now a function of <j> we get: 

Ef nhom {C) = 2kJ dAp± [H±{C,fa) -H ] 2 + 2kJ dA[l - p%\ [H±(C,4>)} 2 , (3) 

where p± are the polymer distributions on the outer and inner monolayers. 

Subtracting Eq. (|^) from Eq. ([!]), and taking into account conservation of polymer, J dAp^ = j dApo, we see that 
the energy gain, AE± hom (C,H ) = E± om (C,H ) - E± hom (C,H ), depends linearly on the spontaneous curvature: 

H ± {CA). (4) 

Thus, for any polymer distribution with J p^H{C,(f))dA < J poH(C, 4>)dA, the inhomogeneity lowers the energy; i.e., 
AE^ lhom is positive and can become arbitrarily large for large values of Hq. The detailed calculation (see below) 
shows that such configurations indeed exist. 

As for the entropy of the system, we assume that the dominant contribution is the entropy of mixing of the polymers 
and lipids. This entropy is larger when the distribution of polymers around the cylindrical bilayer is homogeneous, 
favoring a straight tube. However, it does not depend on the spontaneous curvature. Therefore, if Hq is large enough, 
the energy gain due to AE inhom (C, Hq) is larger than the free energy cost coming from AE^ om and the entropy of 
mixing. In this case, the tube is bent at equilibrium. It remains to be shown that such an equilibrium state can occur 
for reasonable and physical values of the model parameters. For this purpose we now turn to the full calculation. 
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B. Mean Field Calculation 



We neglect fluctuations of the central line curvature C, and correlations between different segments of the tube. A 
section of tube of radius r and length I, with circular cross section having a fixed central line curvature, C, has an 
energy E = E + + E~ , where E + and E~ denote the energy of the outer and inner leaflets, respectively. E ± take the 
form: 

E ± = 2a"Y / H^-4) + ^4] (l G± (<t>i)-4 H °) 2 > ( 5 ) 

where % — 0, . . . , — 1 is the index around the tube and j — 0, , — (1 + Cr cos </>j) is the index along the tube 

(Fig. [9]). afj = 0, 1 are the occupation variables of the polymer on the outer and inner monolayers respectively, and 

G ± ((j>i) = ±.\ ( 1 + ^cTcos^- ) ' ^- e - ( ~' ± ( ( t > i)/ r is the local mean curvature of site (ij). fa — is the angle around 
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the tube as shown in Fig. g. By this convention the curvature of a site on the outer leaflet has the same magnitude 
and opposite sign of the corresponding site on the inner leaflet. 

To find the partition function, as a function of the average polymer occupation, we sum over the polymer degrees 
of freedom: 



{<r±=0,l} 



(6) 



The second term in the exponent is a Lagrange multiplier that allows us to set the average concentrations on the 
membrane, p± = a %l^ ^ to the desired value by adjusting the chemical potentials p ± . N is the total number of 
sites per monolayer. Note that p± is the average of /rr. From this we calculate the free energy: 



F ± (C, p) = -X s Tln(S ± ) + p ± N P± . 
The summation over the polymer occupation degrees of freedom can be carried out, leading to: 



2tt 



d(p(l + Cr cos (j>) In 



(7) 



(8) 



where the approximation ss - d(f> was used. This approximation is valid when the size of a patch is significantly 
smaller than the radius of the tube. We estimate that a/r w 1CU 2 in our system, allowing us to expand Eq. (||) in 
powers of a/r to second order, and evaluate the integral in (0): 
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(l + y±)Vl-(Cr) 2 
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(l + y±)Vl-(Cr) 2 



(9) 



where y ± = exp (/3/x ± - 2f3K'a 2 H§). 

We are interested in finding the free energy as function of the polymer concentration on the membrane, p±, rather 
than the chemical potential p^. Using 



2-nlr 

n * 



we get 



P± = 



dln(S 

d(l3p±, 

2(3ni'a 2 H 2 y ± 



± 



IF 



r(l + y ± ) 2 



O 



(10) 



(11) 



Solving for y^ gives j/^ = p±(l T 2 2 /3/t'ai?o)/(l ~ P±) to second order in a/r. Substituting this into Eq. (^) and 
subtracting the energy for forming the straight tube, F^(0, p±, r), we get the free energy cost of bending a single 
bilayer of radius r within the stack: 



F ± (C,p±,r)-F ± (0,p±,r) 



+4(3 K ' 2 a 2 H 2 p 2 ± ] 



K - 4f3K u a z H£) p± 
1 



VI - (erf 



- 1 



(12) 



The free energy cost of bending the entire tube is then the integral of Eq. (|12|) over the stack. In general, p± 
are functions of r. As we do not know the form of this function, we assume for simplicity that the average polymer 
concentration is the same on all the monolayers of the stack, i.e. p±{r) = po- Variations in p± with r do not 
qualitatively change our conclusions. Under this assumption the free energy cost of bending the tube is 



f(C,po) = 2 I dr[F + (C,p ,r) - F+(0, Po ,r) + F-(C,p 0l r) - F-(O lPo ,r)] = 



2lKtube(po) 



In 



1 + y/l- (Cr y 



(13) 
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where Ktube = -J 1 [4/3k' 2 a 2 i/pp 2 , + (k' — k — A(3 k 12 a 2 H^) p + re] , and d is the bilayer spacing in the stack. 

The logarithm in Eq. (|l3| ) is an increasing function of C. Thus the behavior of the tube is dictated by the sign of 
the effective bending modulus of the tube, Ktube- If it is positive, the minimum of /(C, po) is at C = 0, and tubes are 
straight on average. If, on the other hand, Ktube < 0, tubes form tight coils, since the minimum of the free energy is 
at the maximally possible central line curvature C = 1/t"q. In this case, the free energy of the tube decreases upon 
bending, in agreement with our qualitative argument (see above). 

The typical dependence of K tu be on p for large enough values of H is shown in Fig. |Io| as a solid line. When 
po < pi, Ktube is positive and decreases with po. In this regime the tube is predicted to be straight on the average, but 
with enhanced fluctuations due to the smaller bending modulus. Although we assumed that the presence of polymer 
molecules increases the local bending rigidity of the membrane [k! > k) , their mobility makes it easier to bend the 
tube. 

For pi < po < pi-, Ktube is negative and the tubes form maximally tight coiled structures, pi is therefore a threshold 
occupancy above which straight tubes are unstable. Above P2 straight tubes become stable again. However, this 
regime is probably unreachable in our experiments, since too large a polymer concentration destroys the bilayers. 

For small values of Hq, K tu be is always positive and the model does not predict a coiling instability (dashed line 
in Fig. |l0| ). Therefore, we now check whether the experimental values of the various parameters correspond to a 
regime in which a coiling instability is predicted. Using pipette aspiration pl| we measured the bending modulus of 
a bilayer to be 2k ~ 20 ± hknT. We assume that Hq > lOpm^ 1 because we have observed objects that have radii of 
curvature of the order of, or smaller than optical resolution (~ 0.2pm). We suppose k > k; this is consistent with 
models of composite membranes ]20|-|22]| (although the systems these models describe are different from ours), and 



with experiments (Evans measured k' ~ 2k for membranes with grafted polymers |19|). Putting these estimates into 



Eq. (13) we see that even for small amounts of polymer Ktube can become negative (pi < 0.1), leading to a coiling 
instability. 

The model predicts an inhomogeneous polymer concentration around the tube, which we now wish to calculate. In 
the expression for 2 one can use a (^-dependent chemical potential p^. The distribution of polymer around the bent 
tube is then calculated as follows: 



± dM&«j$}) i (14) 

d(flfl^) 1 + e -2/3^-KG±(0)2+2/3a2 K '(I G ±(0)-ffo) 2 -/3M ± 

The chemical potentials, p/^, corresponding to a particular total concentration, p±, were found numerically by 
plotting the integral of Eq. ([l4]) over 0, as a function of p^ . The value of p^ corresponding to the desired p± was 
then read from the graph. Using it in Eq. ([l4]), we calculated pj. Examples of such distributions are shown in Fig. 

B 

V. IMAGING MODEL 

The model of the previous section predicts maximally curved coils, i.e. C — 1/ro, for pi < po < P2- In order 
to test this prediction we analyzed the experimental pictures in detail. The images obtained from the microscope 
are two-dimensional projections of the viewed object, and include contributions from regions which are out of focus. 
Phase contrast microscopy complicates the interpretation of the images further, since the intensity at a particular 
point is not a monotonic function of the amount of material that the rays of light traverse. This can create apparent 
defects, as illustrated in Fig. [l2], and makes resolving images such as Fig. very difficult. 

In order to test the theoretical predictions, we have used numerical simulations of various coils to derive virtual 
projections, and compared them with the experimental images. The simplest objects we considered are ideal single 
and double helices. We have calculated a geometrical phase diagram, for the different possible types of single and 
double helices, their curvatures and lines of contact, as described in the Appendix. Finally, we have considered 
complex coils reminiscent of the complex structures seen in many of our experiments. 

By generating simulated curves we can see the effects of various imaging parameters on the resulting image. For 
example, when the coil under observation crosses the focal plane, its shape seems to change from a symmetrical 
arrangement to a series of parallel streaks. This is shown in Fig. where an experimental image is compared with 
a simulated ideal helix rotated by ten degrees with respect to the imaging plane. In another case, the image gives the 



impression of a helicity reversal as illustrated in Fig. 12, even though the helix is ideal. A similar effect is seen when 
a myelin figure is slightly curved. This demonstrates that even slices of very simple objects may display complex 
features. Therefore, it is extremely difficult to reconstruct the objects which correspond to the experimental images. 
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Could the tight complex coils observed experimentally (Fig. g|) have maximal central line curvature everywhere? In 
view of the difficulties outlined above, we can only partially answer this question. In order to fully define a curve in 
three-dimensional space, two parameters have to be specified at every point. One choice for such parameters is the 
curvature, C, and the torsion, r. Our theory predicts that C — 1/ro when p\ < p < p%. It does not, however, specify 
what t should be. We investigated various shapes with C w 1/ro and varying r using a numerical simulation. The 
program takes two known functions for C and r, and integrates them to form a three-dimensional curve. Concentric 
cylinders are then drawn around this curve. The resulting object can then be rotated and sliced. Such a slice represents 
to a certain extent the depth of focus of the microscope. The resulting image is then smoothed using a gaussian filter, 
in order to remove the underlying grid, and compared to an experimental image. 

An example of this procedure is illustrated in Fig. [l4|. The best fit that we found to the experimental object shown at 
the top is displayed at the bottom. This is a maximally curved double helix with a complex (though periodic) torsion. 
It may be possible to generate more complex maximally curved coils, similar to those observed (see Fig. ^|e), using a 
non-periodic torsion. However, it is difficult to simulate such objects because of the difficulty in enforcing excluded 
volume constraints. Thus, we can only simulate the simplest structures observed experimentally. The curvatures of 
these structures are consistent with the prediction of the model; i.e., they are maximally curved. 



VI. DISCUSSION 



While our model explains the observed phenomena, other models were proposed in the literature to explain coiling 
in other systems. The discussion is devoted to these other possibilities. 

The coiling of DMPC-cardiolipin myelin figures was attributed to Ca 2+ -mediated membrane-membrane binding 
p5[ . Surface adhesion was also suggested as the dominant mechanism responsible for the coiling of egg- yolk lecithin 
myelin figures . It was proposed that the gain in energy coming from surface adhesion overcomes the elastic energy 
required to bend the tube. Furthermore it was postulated that by twisting back along its length, a rod increases the 
length of the line of contact, thereby increasing the area of contact. An analysis of simple helices shows that this is not 
generally true. This is shown in the Appendix. First, there may be up to four lines of contact, depending on the type 
of double helix. Moreover, the length of the line of contact is not simply related to how tight the structure is. In fact 
the maximum length of the line of contact occurs when the two helices revolve around a large radius, thereby being 
effectively straight. In our system we have found strong evidence that surface adhesion does not play a dominant role 
in determining the morphology. 

A generalization of our model may be relevant to other multicomponent systems. For example, calcium is known 
to induce a lateral phase separation in a lipid mixture containing cardiolipin, thereby creating an inhomogeneous 
distribution. Egg-yolk lecithin is also a mixture of lipids with the same head-group, and a variety of tail-lengths and 
degrees of saturation [Q. 

Another type of curvature model - the Area Difference Elasticity (ADE) model - was discussed in the literature 
p3| . In this approach it is assumed that each monolayer has a certain preferred area. If there is an area difference 
between the monolayers comprising a bilayer, then the bilayer will bend in order to accommodate it. In effect ADE 
and SC both predict the same equilibrium shapes for hollow vesicles. However, ADE cannot explain the coiling of 
myelin figures. This is because bending a tube, so that it is part of a torus, does not change its area. This is true for 
each monolayer in the stack, and in particular, the area difference between layers is unchanged. 

We believe coiling occurs mainly due to the spontaneous curvature induced by the presence of the polymer molecules 
as well as their mobility. Our model emphasizes these aspects, and neglects others such as membrane distortions due 
to the polymer backbone and interactions between polymers. We do not expect these effects to change the qualitative 
behavior of the system. 



VII. ACKNOWLEDGEMENTS 



We thank especially L. Jullien, for his help and encouragement. We also acknowledge useful exchanges with R. 
Lipowsky, E. Moses and S. Safran. This research was supported by The Israel Science Foundation administered by 
the Israel Academy of Sciences and Humanities - Recanati and IDB Group Foundation and The Minerva Foundation. 



7 



APPENDIX A: GEOMETRIC PHASE DIAGRAM 



The requirement that a physical coil does not intersect itself imposes a constraint on the relation between the pitch 
and the radius of a helix. This constraint, as well as the lines of contact, are now calculated for single and double 
helices. Some of the helices studied here were considered also in p^ , p5| . 

1. Single Helix 

The central line of a helical tube can be parameterized as follows: 

Rxit) = (ri cos(i), n sin(i), ht) (Al) 

where n is the radius of the helix around the z-axis, 2irh is the pitch, and t is a parameter running along the curve. 
From this we see that the curvature is C — r\j(r\ + h 2 ), while the torsion is given by t — h / (rf + h 2 ) . 

A myelin figure is then represented by drawing a stack of concentric tubes having a circular cross section around 
this line: 

R(t, 6) = Ri (t) + r \n(t) cos(0) + b{t) cos(0)l (A2) 

where n(t) and b(t) are the normal and binormal to the curve at point t and < r < ro as shown in Fig. ^. The 
demand that a physical tube does not intersect itself yields two conditions: 

1) The curvature of the central line, C < 1/tq. At C = 1/ro two consecutive segments along the tube come into 
contact. From this condition we find that h > y/ r\(r$ — r\). 

2) The pitch of the coil must be such that a segment of the coil does not intersect any previous segment. This 
requirement can be checked numerically by determining the distance between any two points along the curve, and 
demanding this distance to be greater or equal to 2r$. The result is another constraint on h(j*o, r\). Helices for which 
this distance equals 2tq have a contact line which is a helix with the same radius and pitch as the central line, but 
shifted by r along the line of symmetry. The allowed parameters are summarized in Fig. ^|. 

2. Double Helices 

We consider a double helix as two single helices, having radii r\ and r 2 , that have at least one continuous line of 
contact. If there is no such contact then they are simply two unconnected single helices. For simplicity we assume 
that both tubes have the same radius ro around their respective central lines, and that these central lines wind around 
the z-axis. The line of contact is itself a helix with the same pitch, and with a radius r 2 < r con < r%. The constraint 
of having a distance of 2r$ between the two central lines allows us to find r\ as a function of r 2 and h. 

There are two types of double helices. Type I double helices have h > y 'ri(2r^ — r^J. In this case, the two helices 
wind around their line of contact and their radii obey t\ = 2rg — ri. Double helices of type II obey the condition 
h < y ri(2rQ — ri). In this case, the coils are intermingled and wind around a hole in the middle. Type II double 
helices have at least two lines of contact (an example is shown in structure d of Fig. [H]. This is summarized as a 
phase diagram in Fig. [L6l 



We now find the lines of contact between the two coils of a double helix. Using Eq. ( Al) we define the central lines 
of the two strands: 

-Ri(^i) = (n cos(£i),ri sin(ii),/iii) 

^2(^2) = (^"2 coa(t 2 + 77), r 2 sin(i 2 + 77), ht 2 ) . (A3) 

The phase r\ is r\ = n for double helices of types I and II, and for loose helices 77 is determined by the requirement 
that there exists at least one line of contact. The line of contact is then given by D(ti,ti) ~ Ri(ti) — Ri(ti). From 
the requirement that D _L dRi/dti and that \D\ = 2ro we find that the line of contact is a helix with the same pitch 
as the other two helices, and a radius: 

1 



2nr 2 cos(A), (A4) 
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where A = Jb/d - 2 ± 2y/l - b/d + 1/d 2 , d = h 2 /{nr 2 ) and b = (4rg - r? - r|)/(r 1 r 2 ). 

Double helices of type I have one line of contact, while type II helices typically have two lines of contact. At the 
dotted line in Fig. |l^ the two lines merge. Type II helices lying on the line of self intersection also have lines of self 
contact, for a total of four lines of contact. The maximal length of contact lines occurs when one helix lies within 
the other and their radii go to infinity. In this case the total length of contact is four times the length of the central 
line of each helix. Thus, formation of a tighter helix does not necessarily mean that the line of contact increases, not 
to mention the area of contact. The area of contact between two cylinders depends on the angle at which they come 
into contact. Moreover, when there are two lines of contact close together (for example near the dotted line in Fig. 
|l6|) the areas of contact must overlap to some degree. 
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FIG. 1. Image of a myelin figure at very low polymer concentration. The scalebar represents 20 /im. 



FIG. 2. When the instability begins at the tip, the tube loops back upon itself forming a globular structure. Such a structure 
is seen here (a) (b) 83 (c) 109 (d) 155 (e) 221 (f) 371 seconds after onset of bending. The scalebar represents 10 fim. 



FIG. 3. A sequence of images depicting a single helix being mechanically stretched, and returning to a maximally curved 
configuration. The helix behaves like a spring, responding to the stretching force by elongating uniformly. When the force 
is removed the coil retracts. This behavior suggests a restoring force, rather than surface adhesion. Times are: (a) before 
stretching, and (b) 67 (c) 102 (d) 159 (e) 215 (f) 325 (g) 393 seconds after initiation of stretching. Snapshots (d-g) were taken 
after the force was removed. The scalebar represents 10 ^m. 



FIG. 4. Fluorescence image of a double helix. The scalebar represents 10 /im. 



FIG. 5. Formation of a complex coil. The tube becomes unstable locally, forming a hairpin which gradually curls up. 
The intervals between snapshots (a-d) are 45 seconds long. Structure (e) was observed 13 minutes after onset. The scalebar 
represents 10 /im. 



FIG. 6. A double helix coming into contact with another myelin figure. As the coil grows, it pushes the other figure aside. 
Fifteen minutes after initial contact, the two myelin figures are no longer touching. This is another demonstration that adhesion 
is not important in this system. The scalebar represents 10 /im. 



FIG. 7. Fluorescence image of a myelin figure showing that there is polymer inside. Inset: fluorescence intensity along the 
section. The scalebar represents 10 /im. 



FIG. 8. Schematic representation of the heuristic argument. Starting from a straight tube with the same polymer concen- 
tration on the inner and outer monolayers (a), the cost of bending the tube keeping p(<j>) homogeneous (b) is E^ om - Allowing 
the polymer to diffuse around the tube (c) to the configuration shown in Fig. [ll] lowers the energy by A_E* hom , but also lowers 
the entropy. The black (white) color in (c) corresponds to areas of small (large) polymer concentration. Both Ehom and the 
entropy are independent of Ho, while Er nhom is linear in Ho- Thus for large values of Ho it is preferable to bend the tube. 



FIG. 9. A small section of a coiled stack of concentric membranes showing the directions of the the normal, n, and the 
binomial, b. Each monolayer, having a radius r is divided into patches with running indices i and j. 



FIG. 10. The effective bending modulus of the tube, K tu be, is parabolic in the average occupancy, p. We have used the 
following values of the parameters: k = lOfcsT, k — 2k and ro = 5/tm. We find that Ktube depends on a and Ho only through 
the product aHo- The solid curve represents Kfube for aHo = 0.3. When aHo is large enough (aHo > 0.19 for the values of k, 
and k' we have used), Ktube < between p\ and pi. For smaller values of aHo, K-tube > for all values of p. The dashed curve 
corresponds to aHo = 0.16. 
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FIG. 11. The local polymer concentration as a function of the angle around the tube is shown for two average polymer 
concentrations: p — 0.5 (solid line) and p — 0.1 (dashed line). The following parameters have been used: r/a = 100, rC = 0.99, 
aHo = 0.3, k — WksT, k — 2k, and p = 0.5 (solid line), p = 0.1 (dashed line). 



FIG. 12. Top: Fluorescence image of a double helix showing that the helicity seems to reverse where the coil intersects the 
focal plane. This is an artifact of the imaging geometry, as can be seen for the theoretical ideal double helix tilted at 15° at 
the bottom. The scalebar represents 10 pm. 



FIG. 13. A comparison between an experimental helix (left) with a model ideal double helix (right). This image shows the 
effect of tilting the coil by 10° with respect to the imaging plane. Notice the cut has a symmetric appearance at the bottom, 
where the imaging plane slices through the center of the coil, shifting gradually towards an array of parallel smudges at the 
top. The scalebar represents 10 pm. 



FIG. 14. Top: experimental image of a double helix. Bottom: cross section of a theoretical double helix having maximal 
curvature (C = l/^o), and sinusoidal torsion, r = 2 + 2sin {ff [i+ 1.8sin(^j 2 - + 4^)] }. The central line of the image is shifted 
by ro below the imaging plane. The scalebar represents 10 pm. 



FIG. 15. The line of self-contact for single-helices, ro is the outer radius of the tube. 2nh and r\ are the pitch and the radius 
of the central line of the helix, respectively. The requirement that a physical coil does not intersect itself gives two conditions. 
The dashed line corresponds to the first condition, while the solid lines corresponds to the second condition, as explained in 
the text. 



FIG. 16. Geometrical phase diagram for double helices, ro is the outer radius of the tubes and 2irh is the pitch of the central 
line of the two helices, n and r2 are the radii of the central lines of the two helices. Our convention is that T2 < r\. There are 
two types of geometries, separated by the dashed line in the figure. Type I: two helices winding around each other and type 
II: interlaced helices revolving around a central hole. The lower bound is the limit below which a helix intersects itself. The 
top-right bound is the line of symmetry, where both helices have the same radius. 



FIG. 17. Simulations of representative ideal theoretical helices. Letters are keyed to the diagram in Fig. [16|. Most of the 
double helices observed experimentaly resemble type (c), although some of them are more symmetric than the example shown. 
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